Estimating total momentum at finite distances 
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We study the difficuhies associated with the evaluation of the total Bondi momentum at finite 
distances around the central source of a general (asymptotically flat) spacetime. Since the total 
momentum is only rigorously defined at future null infinity, both finite distance and gauge effects 
must be taken into account for a correct computation of this quantity. Our discussion is applica- 
ble in general contexts but is particularly relevant in numerically constructed spacetimes for both 
extracting important physical information and assessing the accuracy of additional quantities. 

PACS numbers: 04.25.Dm, 04.30.Db 



I. INTRODUCTION 

The possibility of studying highly dynamical, compact 
systems by the observation of gravitational waves with 
the new generation of gravitational wave detectors is 
about to become a reality. Efforts at the experimental 
front have allowed the coming on-line of several highly 
sensitive detectors which will be followed by even more 
refined onesP, [2, Confronting the output with theo- 
retical models will shed light to a number of spectacular 
events which have so far remained out of reach. In fact, 
the analysis of the detectors' output has already provided 
upper bounds to several systems, giving insight into the 
possible source or location of a gamma ray burst event 
and constraining important physical parameters in the 
crab pulsar[5j. A detection and subsequent analysis will 
go further than this by providing key information of the 
system that produced it, if accurate models are available 
to tie theory with observation. For binary black hole 
systems, concentrated efforts are now focusing on provid- 
ing effective templates covering from early orbiting stages 
through the merger andpost-merger stages and assessing 
their accuracy [a^S M, P|- Intimately related to this lat- 
ter issue is the suitable identification of physical param- 
eters at each of the different stages (treated by different 
techniques) and the extraction of sought after informa- 
tion. When studying an isolated system in general rela- 
tivity, there are several quantities that have strict physi- 
cal meaning only in the asymptotic region; for example: 
gravitational radiation fields, total momentum, total an- 
gular momentum, etc.; which are defined at null infinity. 
Numerical approximations however, often are not able to 
reach the asymptotic region and so the relevant quantities 
must either be computed at finite distances or propagated 
somehow to infinity for its unambiguous computation. 
In a previous article [l0|. we have analyzed the difficul- 
ties that appear when estimating gravitational radiation 
at finite distances and how to remove them. In this arti- 
cle, we discuss the subject of estimating total momentum 
also at finite distances from the sources. This is useful 



to extract further physical information and to check the 
numerical solution itself. In particular we concentrate 
on the total Bondi momentum which is an important 
physical quantity. For example, it is directly related to 
the computation of velocity recoils produced by kicks in 
collisions of compact body systems, or in the gravita- 
tional radiation recoils produced in anisotropic collapse 
of stars, etcfH IH, El, Hi, M, MMMMM- The 
prevalent approach to compute these quantities within 
a simulation is to do so through a temporal integration 
of the Bondi momentum flux. The Bondi momentum 
on the other hand needs only be evaluated at the initial 
and final times for this purpose. Furthermore a balance 
law relates the difference between the Bondi momentum 
with the integral of the momentum flux, which can be 
employed as a further check of the implementation. 

In what follows we review the Bondi momentum cal- 
culation at infinity. This will serve as a motivation for 
the rest of the work, by addressing the limitations that 
finite-distance calculations must address to reduce the 
associated errors when attempting such a task. The no- 
tion of total momentum is normally associated to an 
integration on an asymptotic sphere at future null in- 
finity of an asymptotically flat spacetime. The integra- 
tion employs the so called four-translations of the Bondi- 
Metzner-Sachs group (BMS), and an appropriate inte- 
grand. A typical expression in terms of standard notation 
is 
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r(«'° + cr"^°)d52; 
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where P = (1, sin(0) cos(0), sin(0) sin((/>), cos(0)) when 
expressed in standard angular coordinates {9, 0) and all 
these quantities are calculated in terms of a Bondi coordi- 
nate and tetrad systemfs^j- In particular, recall that the 
Bondi mass M is the timelike component of this vector; 
namely 



(2) 
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The Bondi momentum can also be expressed in terms of 
the Psi supermomentumfiH = ^2 + 0"°^° + by 

r^-dS^. (3) 
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This supermomentum has a couple of useful properties; 
namely, it is real, ^' = 'J, and its time derivative is simply 



(7 (T . 



(4) 



Thus, one can see that the time variation of the Bondi 
momentum, in terms of time Bondi coordinate, is just 
given by 
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As we discussed in [13], unless further structure is 
added to 3 + 1 implementations, these can not access fu- 
ture null infinity to calculate useful quantities and must 
resort to do so at finite distances. This might not only 
bring about finite size effects but also difficulties in adopt- 
ing the appropriate frame -analogous to a Bondi one- 
since the advantageous coordinates at the evolution level 
need not agree with Bondi-type ones. In \ld \ we exam- 
ined how these issues influence the calculation of the ra- 
diation and a way to address them. We next carry out 
a similar analysis for the momentum. Certainly, as the 
underlying issues are the same, we rely on our discussion 
and refer the reader to that work for further details. For 
the sake of completeness however, we include here briefly 
the list of basic assumptions required in the calculation of 
relevant quantities at future null infinity (further details 
are presented in [l^) : 

• Peeling is assumed. 

• Outgoing null hypersurfaces, parameterized by u 
intersect T'^ (future null infinity) defining a se- 
quence of S*^ surfaces. 

• Each of these surfaces is conformal to a unit sphere 
metric; off this surface (into the spacetime) the de- 
parture from it is of lower order. Namely the an- 
gular metric in a neighborhood of 2^ can be ex- 
pressed, as gAB = r'^hAB = r'^iqAs/V'^ + cab/t + 
0{r~^)); with V a conformal factor, qab a unit 
sphere metric, and r a suitably defined radial dis- 
tance. 

• A null-tetrad rt°, m", m°} satisfying i°-na = 
—171°" fha = 1 (with all other products being 0). 

• Using standard ^,'23'] conventions for the Riemann 
tensor and spinor dyad, particularly useful scalars 



obtained from them are. 
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A suitable (Bondi type) expansion in terms of 1/r 
with coordinates chosen such that V — 1, g^^ ~ 
1 and g^ji^ — [x^ labeling angular coordinates 
at u = const, r — > 00) gives rise, in particular, 

to 



9*^ 
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where the supra-index indicates leading order 
in an expansion in the radial coordinate and 9 is 
the edth operator[22] of the unit sphere. However, 
if one does not adopt a Bondi frame, the previous 
equations become much more complicated. 

Working with finite size regions one would like to esti- 
mate the asymptotic fields in terms of null tetrads based 
on a choice of coordinate system. The nature of the 
numerical work will suggest the choice of coordinates 
-adopted to simplify the simulation- and consequently 
some natural tetrads that can be defined. These null 
tetrads however, will not be in general Bondi tetrads. 
In order to make these estimations, an extraction world- 
tube is assumed to be far enough to calculate the asymp- 
totic quantities. See figure [T] 




FIG. 1: . The extraction world-tube F and the associated 
null tetrad and coordinates. The nearly vertical curves de- 
note the time evolution of the angular coordinates on the 
world-tube. The tetrad null vector i points outwards, while 
n points inwards. Complex null vectors (m, m) are tangent 
to the spheres contained in the world-tube. 

The metric is expressed in terms of a null tetrad 
through 
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(13) 



The null tetrad can be related to a dyad of spinors 
(o^, (.'^). The relation among the null tetrad vectors and 
a spinor dyad is given by ^" <^=^ 0^0^ , m" ■^=^ o^l^ , 
771° i^o"^' and n°- ^=> l^l"^' . 

Therefore, the transformation of the tetrad frame 
can be obtained, as we do here, from the transfor- 
mation of the dyad frame[3lj. We will also make 
use of a null polar coordinate system {x^ , , , x^) = 
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(m, r, (C + C), i(C — 0); where r has the meaning of a ra- 
dial coordinate and C) angular coordinates. 

Using the numerical time coordinate i, we construct 
the coordinate system {u, r, Q in the following way. We 
define F to be the timelike surface given by (r = i? = 
const.). On F, a null tetrad {€, n,TO,m} can be defined 
in the following way. Let the null function u be such that 
on F one has u = t; and £ — du everywhere. The function 
u is chosen such that the future directed vector £ points 
outwards with respect to the 2-surfaces {t — const., r = 
R); which are, topologically, two dimensional spheres. 

Then, the complex vectors m and m are defined to be 
tangent to the spheres {u = const., r = const.). Further- 
more, one can choose m to be proportional to and 

TO to be proportional to in the asymptotic region for 
large r. The remaining null vector n° is set by requiring 
the standard normalization conditions. This is the way 
in which the coordinate system (it, r, (, is constructed. 

Now, in order to compute the total momentum one 
must be able to express the different Bondi quantities 
appearing in eq.([T]), in terms of analogous quantities con- 
structed from the tetrad {^°, n", m", m"'} and coordinate 
system (m, r, In the next sections we discuss how to 
do so. 



II. BASIC STRUCTURE OF 
ASYMPTOTICALLY FLAT SPACETIMES 

A. Asymptotic behavior of the tetrad's 
components 

Let us use the coordinate system {u,r,(,C)- Keeping 
(C = const.) and {( = const.) on the null hypersurface 
(u — const.), one can see that incrementing r one moves 
along a null direction. Since £ is contained on the hyper- 
surface {u = const.), one deduces that £ is proportional 
to We wiU prove below that 



A natural null tetrad is then completed with. 
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dr 



with i = 2, 3 and components 5*, U and X* are: 

t2 / 1 \ f3 
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with 



(20) 



where V = V{uX,0 and the square of Pq — ^^-^7^ is 
the conformal factor of the unit sphere; 
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U = rUoo + Uo + — + o(\ 



(21) 



and the other components of the vector have the 
asymptotic form 

^2 = Xl + ( i j , X"^ ^Xl + O ( i j . (22) 



B. Coordinate and tetrad transformations 

With the objective of computing the total momentum, 
one must transform to an asymptotic Bondi coordinate 
(m, f, C, C) and tetrad frame {£, n, m, m). 

The transformation is of the form 



u = a(M,C,C) + 



"i(",C,C) 



+ O (r") 



C = Co(",C,C) + o 



(23) 
(24) 

(25) 



with d > 0. 

Let us note that in a Bondi system, the angular co- 
ordinates (C, C) can be chosen as the stereographic co- 
ordinates of asymptotic spheres. If one further assumes 
that C is a stereographic coordinate of the spheres [t = 
const., r = const.); which are conformally related to the 
Bondi coordinates, it is only necessary to consider an 
angular transformation of the form[To| 



C = Co(^,C) + o 



(26) 



The contravariant metric components for the stan- 
dard Bondi like coordinate system are given by equa- 
tions (3.13-18) of [1^; whose inverse is given by equations 
(3.19-24) of the same reference. 

For the general tetrad the only difference is that now 



9ur 

has not necessarily a unit value. 



(27) 
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In order to deduce the relations among the dyads let 
us start from the null tetrad defined by p^ . 



(D 



dr 



(28) 



together with eqns. p^ - p^ . Then the inverse metric 
is given by 





= 0, 


(29) 




= A, 


(30) 




= 0, 


(31) 




= 2AU, 
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While the metric is given by 
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— A 

grr = 0, 
gri = 0, 
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(35) 

(36) 
(37) 

(38) 

gui = -gi]X^, (39) 
g^J = ig'')-' = -de.kejiie^' + e'CO; (40) 

with ij, fc, Z = 2, 3, d = det(gij), = -e^j and £23 = 1- 
In particular, defining the quantity 



gnu = -2^ + X'X^gij, 



one has that 



d = 



(41) 



(42) 



The coordinate transformation induces a tetrad trans- 
formation as described in [T^j ■ The main resulting equa- 
tions are: 



•1= ^± 
AO' 



(43) 



which determines a relation among w, a and (7"^, where 
a superscript 'O' means leading behavior for large r. 
The function ( must be chosen [To| such that 



C + x°^ Q = 



(44) 



so that angular coordinates remain constant along the 
generators of future null infinity. In particular if to lead- 
ing order gui = (X* = 0), one could adopt initial con- 
ditions so that C C- III the general case however, the 
coordinates will "shift" around the world tube and the 
explicit transformation must be taken into account as 
described in [26( |. 



The transformation of the vector n is given by 
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(45) 



Similarly, the asymptotic behavior of the vector m is 
given by 
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where is the vector n evaluated at future null infinity; 
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(47) 



and therefore an operator acting on functions. 

Since the angular part of the metric expressed in terms 
of the new null tetrad must coincide with the angular 
part of the metric expressed in terms of the original null 
tetrad, it follows that 



^0 ^ C C = 1 



So that 
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and therefore the factor of m is 
wPoV 



CcCc ; 



(48) 



(49) 
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(50) 



Then, in particular, for a Bondi system one has 



where we are using that the tilde system is of Bondi type. 
Equation (|5ip determines a. 

From the transformation expressed in equations (j45p 
and (|46p and the requirement that the spinor metric sab 
be invariant, one can deduce the transformation for the 
spinor dyad 




(52) 
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and 
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(53) 



by 



The regular dyad at future null infinity is then given 




(54) 



and 



5^ 




where 



B 



and 
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(55) 



(56) 



(57) 



Notice that one could generalize (|52|) and (|53p to in- 
clude higher order terms in an asymptotic expansion of 
the form 



1 
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(58) 



and 



B(o^(l + (5) + fiC(l+7)i^); (59) 



however its contribution to the (spinorial) metric is of 
higher order as shown by the following expression 

= [(l + (5)-r!/iC(l + 7)]e'4^; (60) 

which in turn implies, 

<5 = 0/iC(l + 7). (61) 



III. CALCULATION OF THE BONDI 
MOMENTUM 

In order to compute the Bondi momentum, as ex- 
pressed in eq. HI) , one needs to evaluate the Bondi quan- 
tities s-^d ^''^ ™ terms of the quantities obtained 
in the more readily accessible quantities to numerical im- 
plementations. 



A. Transformation of 

We can now easily calculate the component ^'2 of the 
Weyl tensor, in leading orders, with respect to the new 
null tetrad, obtaining 

~A~B~C~D 



B^ 



*2 - 2 
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(62) 



Note that the asymptotic transformation of the dyad 
{o^,i^} can be re-interpreted as a combination of the 
null rotations of types II (rotation around n"-) and III 
(boost / spin) in the standard Newman Penrose formalism 
(see Appendix A). If one makes use of these rotations to 
compute the transformation of the complete scalar '^2 an 
expression similar to eq. (j62p is obtained but without the 

factor ((7°^ ■ This factor is present in eq. ([5^ since 

we are computing the leading order behaviors, which due 
to the assumed peeling property obey, 



*2 ^ 1 



*2 
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Then, taking into account the relation 



r + O (r°) 



(63) 



(64) 



(65) 



if follows that the aforementioned factor will be present 
in the transformation of 



B. Transformation of a° and cr° 



Let us recall that the shear is defined by 



a = 0^1^' o^Vaa'Ob = ^ + O ( ^ ) . (66) 



Only the leading order transformation of the dyad is 
needed to calculate the transformation of the leading or- 
der behavior of the shear. 
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Then, one has 



IV. EXAMPLES 



{wBf 1^^^ [a-n rn{nc) ~ ncn{nc) 

ncr + 2nc {13 - nC(') 

+ n'^c^ i-p' - ncK') + o{n^)] ; 



(67) 



where k — 0. 

The details of the calculations leading to ctq are given 
in the Appendix B. One finds 



ao~dvC + 2^ QyA - C n°(C) 



(68) 

and where C is recognized as a quantity of boost-spin 
weight {p,q} type {2,0} [13]. 

In the calculation of the Bondi momentum, one also 
needs the time derivative of the shear that is given by 



CTO 



du 
1 



n"(a) 



Sc. 



(69) 



(TO - SyC + 2^ dvA - C'n°(C) 



Then, we see that by expressing the Bondi scalars 5'2, 
do and ctq in terms of ^t^, <7o and (tq, we get non-trivial 
factors and terms containing ^tg, ^'4, C, a, V and 
some of their derivatives. This is an analogous situation 
as that observed in [l^l for the radiation field '^^ where 
these correcting terms will have in general non-trivial an- 
gular dependence which would affect relevant quantities. 
These would be important, for example, in computations 
of kicks observed in collisions of different masses or spin- 
ning compact objects where the merged black hole moves 
through the computational grid and, consequently, the 
extraction sphere will be non-centered introducing extra 
multipolar structure. 

For completeness of the discussion we point out how 
to calculate the auxiliary functions: the conformal func- 
tion V, was explained in our previous article ilOj, the 
angular transformation ( is determined by the angular 
components of the vector n, while a is determined from 
the knowledge of gur, C ^• 



In what follows we discuss two examples to illustrate 
the effects brought by non-adapted tetrads in the calcula- 
tion of the momentum. The first example shows how the 
computation can be affected already at future null infin- 
ity and without even introducing a 3-1-1 description of the 
spacetime. The second example illustrates how coordi- 
nates not adapted to a Bondi system in a 3+1 description 
can give rise to results with spurious gauge dependence. 



A. Robinson- Trautman metrics 

As a first example of applications of our expression 
for the Bondi momentum, let us study the Robinson- 
Trautman metrics. These metrics are vacuum solutions 
of Einstein's equations containing a congruence of diverg- 
ing null geodesies, with vanishing shear and twist [27|. 
They can be written as[28j 



-2Hr + K 



M{u) 



du +2 du dr 



jdC dC, 



(70) 

where P = P{u,CX), H = ^, K = AlnP, a doted 
quantity denotes its time derivative and A is the two- 
dimensional Laplacian for the two-surfaces u ^constant, 
r =constant with line element 

dS^ = 4? rfC dC 
P"^ 

We can describe this line element in terms of the line 
element of the unit sphere; this is done by expressing P 
in terms of P = V^(u, C, C)^o(Cj C)i with Pq the value of P 
for the unit sphere. The function V^(u, C, C), must satisfy 
a fourth order parabolic differential equation 



(71) 



where we fixed the freedom of redefining u in such a way 
that M{u) is actually a constant. This equation is known 
as the Robinson- Trautman (RT) equation. If / denotes 
the vector field that generates the null congruence, then 
la = du, and P = (^) . By completing the tetrad in 
the usual way with the three null vectors, 



d_ 

du 

V2P 

r 

V2P 



rr K M\ 

Hr \ 

2 r 



d_ 
d_ 



8-r) ' ^''^ 
(73) 

(74) 



we have a frame adapted to the geometry. We call this 
coordinate/tetrad system RT frame. 

The relation between the RT coordinate system and a 
Bondi system {u, f, C, will be of the form given by eqs. 
(IMD and (I2S1). 
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Now, let us compute the Bondi linear momentum 

P" = / + a°^] r dS^; (75) 

47r J \ " ou J 

in a given u ^ Uq = const, RT section, and let us 
write this P° in terms of an RT frame. The relation 
between Robinson- Trautman and Bondi quantities are 
given by eqs. ([62|l . ([68|l and ([69]) . But in order to sim- 
plify the expressions we will use the freedom of BMS 
group in such a way that the "origin" of the new Bondi 
system u = coincides with the RT section uq, where 
we wish to compute the momentum. That is, the re- 
lation between the Bondi and RT coordinates is such 
that a(uo,CjC) = 0- This implies in particular, that 
3ya(uo, C, C) = ^va{uo, C, C) = 0. 

For Robinson- Trautman metrics we have: 



\I>o 

^2RT 
9ur 

a 



-M, 



a 

1, 
V. 



(76) 
(77) 
(78) 
(79) 



then in terms of these quantities, and using the eqs. ([62|) . 
([68|) and ([69|) . we have that the Bondi momentum reads, 



pa 
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(80) 



Here it is important to emphasize that the absence of the 
correcting factor would imply an erroneous constant 
mass result. 

Note additionally, that if we wish to compute the Bondi 
momentum at two different RT times using eg. ([SO)) , we 
will obtain two Bondi momenta calculated with respect 
to two different Bondi systems which are related by some 
BMS transformation. In order to compare the second 
momentum with the first, one should express it with re- 
spect to the first Bondi system. 



B. Vaidya Solution 

As a second example let us study the Vaidya metric 
which represents a spherically symmetric solution radiat- 
ing a null fluid. The usual representation of this metric 
is 



ds^ 



M{u) 



du^ + 2dudr - r'^dfl'^ 



where 



dn^ = d0^ + sin^ 



(81) 



(82) 



and {u, r, 6, cf)} define a null polar coordinate system. The 
mass of this spacetime at a given time u is given by M{u) 
which decreases with time. Notice that the coordinate 



system used to express the line element is naturally of 
Bondi type with an associated tetrad given by 



1 



V2f 
1 



e 
9? 



sitlO 
i 



-5? 



V2f \ " sin 6* 
Thus leading order behavior of is 



(83) 
(84) 

(85) 

(86) 



(87) 



in terms of an asymptotic regular frame at future null 
infinity. 

We can now compute some geometrical quantities per- 
tinent to the present discussion with this tetrad. In par- 
ticular we obtain for the Weyl's scalars and the shear 
are 



4^2 



■03 = 1p4 

M{u) 



0, 



(88) 
(89) 



thus, in particular, the only non-vanishing component 
of the leading order behavior of Weyl's scalar is -02 — 
—Mill), and the null congruence determined by P is a 
free shear congruence (ct = 0). If we compute the Bondi 
momentum to this metric we obtain the well know result 

P-=-^ J{i^"2+^"~hi-dS^^{M{u),0,0,0y, (90) 

where a is a numerical index. This relation shows that 
M{u) is the Bondi mass, and that this particular Bondi 
system is at rest with the source (the spatial part of the 
Bondi momentum is zero). 

So far the discussion has been centered around a null 
coordinate system. However, most numerical simulations 
adopt a natural 3 -t- 1 coordinate system and it is with 
respect to this system that relevant quantities are cal- 
culated. In order to make contact with such task, we 
must re-express the line element in terms of a 3 -I- 1 
coordinate system through some given transformation 
{u,f,9,(j)) (r, i?, 9,$). Naturally different transfor- 
mations can be envisaged and in particular ones that 
will be consistent with a Bondi system. However, our 
goal here is to illustrate the issues discussed in the pre- 
vious sections in a simple scenario where the coordinates 
do not necessarily conform to Bondi ones. We therefore 
introduce the following transformation. 
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T = A({t)- 

r 

R 



A(i)(w)' 



A(i)(u)' 
6 = 0, 
$ = 



(91) 

(92) 

(93) 
(94) 



where A^-"^^ = ^ and in general we will have A*^") = ^5^- 
Incidentally, let us note that A is playing the role of 
the inverse of the transformation a presented previously. 
Therefore we have. 



dii = 
df — 



1 



-AdT-dR) , 



A(i) 
A(2)i? 



(95) 



A(i) 

de ^ dQ , 

dd) = d<^ . 



+ ( A(i) - 4^ Mi? (96) 



A(i) 



(97) 
(98) 



With this transformation, wc can construct a "stan- 
dard" null tetrad in the usual way as in numerical ef- 
forts, i.e., from a unit vector A'^" orthogonal to the sur- 
face T = const., and three vectors cJq, 9|} obtained 
on the sphere defined by the T = const., R — const, by 
Gram-Schmidt procedure. The result is. 



1 



l"" = —{N^' + R") 



V2 
1 

V2 
1 



V2N 



[dr + Or] , (99) 



— {N" - R") 

[dr ~ {2N^ - l)dB] , 



V2N 
1 



V2A(i)i? 
1 

V2A(i)i? 



sin 9 

i 

sin 9 



9* 



(100) 

(101) 
(102) 



where 

N= J2 



(A(i))2 



Miu) 2A(2) \ , ^ 



and where r and u must be thought of as functions of 
{T,R\. If we compute relevant scalar quantities with 
this tetrad, the only non vanishing scalar is '4'2 which is. 



M{u) _ 
(A(i))3' 



(104) 



where one finds an additional non-trivial factor which 
is completely coordinate dependent. In order to obtain 
the correct expression we need to take into account the 
corrective factors that were discussed in section Hill For 



this, assuming one starts with the 3-1-1 version of the 
Vaidya metric induced by the coordinates (T, i?, 6, <l>), 
and following our discussion we need first a null coor- 
dinate/tetrad system naturally induced by these coordi- 
nates, namely 



u = T-R, 
r = R, 

e = e, 



(105) 
(106) 
(107) 
(108) 



We can now construct the 'naturally' induced tetrad 
{Z", n", m", to"} to leading order and compare this with 
the Bondi tetrad {Z°, n°, m°, to°}: 



la = idu)a, = 9,-? - AW^'' 



1\ 1 



(109) 



= d^^ + O^-j^^n^ + O^-j, (110) 
" V^ri'^^^o'^)-^'^ (111) 



TO 



Then, from eqs. (ll09p and (|110|) . we see that 

gur = 1, a = -^^y. 
If we compute ip2 with this tetrad we obtain, 

V2-^2 - (^(1))3- 

Then, employing equation (I62p . we get 



(113) 



(114) 



(115) 



which yields the correct result. 

As a final comment, let us note that ^2 = V^2- This is 
due to the fact that the relation between the two tetrads 
is 



(116) 























fh'" 





(117) 
(118) 
(119) 



and therefore the ^/2N factors cancel. 



V. FINAL COMMENTS 

The calculation of total momentum at finite distances 
can be employed to extract valuable information from a 
numerically obtained spacetime both to obtain relevant 
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physical quantities and check related results obtained in 
the simulation. In particular, the final black hole re- 
coil velocity can be readily obtained and the total energy 
radiated during the length of the simulation. Standard 
calculations of kicks in numerical relativity have relied 
on integrating the flux of energy over a period of inter- 
est; alternatively, it could also be computed simply by 
the difference between initial and final space-time mo- 
menta; provided that gauge issues are addressed so that 
they do not negatively influence either the calculation or 
the comparison between the results of the two methods. 
Our work provides a way of dealing with possible prob- 
lems and sets up a common framework, consistent with 
Bondi's construction, to remove ambiguities from the cal- 
culation and facilitate the comparison of results obtained 
by different numerical implementations. 

It is probably worthwhile to note that we have been 
concerned with the notion of total momentum at future 
null infinity, as one is interested in the relation between 
total momentum and outgoing radiation which can carry 
momentum. For this reason our approach has been to 
imagine that topological closed 2-surfaces at finite but 
large distances can be understood as limiting spheres in 
the vicinity of future null infinity. However, once one has 
a finite closed 2-surface, one can wonder whether to re- 
gard it as a limiting sphere in a vicinity of past null infin- 
ity instead. With this point of view, one should relate the 
difference in the total momentum at two different finite 
spheres with the incoming radiation; which in this case 
will be characterized by the leading order behavior, at 
constant advanced time, of 5*0. Consequently the notion 
of total momentum at finite spheres will give different 
quantities whether one regards it as limiting spheres to 
future null infinity or to past null infinity. This brings the 
issue of whether one should consider both contributions 
when dealing with finite spheres. The answer is nega- 
tive if one wants to describe what is observable at very 
long distances from the sources; since, for all practical 
purposes, an observer will be situated at future null in- 
finity. Therefore, the notion of total momentum and its 
flux should be those naturally appearing at future null 
infinity. 

To be more concrete on this point, let us consider a 
sphere 5 at a finite distance belonging to a timelike hy- 
persurface F; and which can be regarded as one of a se- 
quence of limiting spheres to future null infinity. Then 
in a numerical calculation one constructs these spheres 
far enough and aims to have enough control on bound- 
ary conditions so that incoming radiation through F is 
negligible. In this scenario, the estimation of total mo- 
mentum by using the limiting behavior of fields in the 
interior, gives accurate information of the total momen- 
tum of a limiting sphere at future null infinity after gauge 
issues are appropriately handled. 

Last, we find it worthwhile to stress that the calcula- 
tion of momentum changes from the asymptotic fields at 
future null infinity does take into account all contribu- 
tions in the interior of the spacetime, including incoming 



radiation that ultimately might cross the horizon of a fi- 
nal black hole. The role of this radiation in the computed 
kick velocity is accounted by the measured variation of 
spatial momentum and black hole mass which changes 
as energy flows into the horizon. Thus, the expression 



kick 



APyP° accounts for all relevant contributions. 
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APPENDIX A 

The transformation between the dyad {o^,t^} and 
|o'^,i"^} induces a transformation on the corresponding 
null tetrads. If we work with regular quantities at future 
null infinity we have. 



m 



(AI) 
(A2) 
(A3) 



and similar expression for its 'tilde' version. Using the 
relations ([5^ and and having into account that the 
conformal factors are related as 



(A4) 



we find 



— i^r + Cm'' + Cm" + CCn'^ ) , (A5) 



= — (m° + Cn°), 
w 



(A6) 
(A7) 



where in order to simplify the notation we have intro- 
duced 



-I, 



(A8) 



From this transformation of the regular null tetrad we 
can see that it is originated as a combination of transfor- 
mations of types II and III. The factor nP{a) represents 
locally the action of a boost, similarly, C the action of 
a null rotation around n and e*^ a spatial rotation. The 
factor w comes as a result of both null tetrads belonging 
to the same conformal metric class (see eq. (jA4p ). 
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APPENDIX B 

The calculation of cto is conveniently done referring to 
a 'starred' tetrad * which is defined in the following way. 
Let us first define a * coordinate system such that 



d 
dr* ' 



(Bl) 



in other words, r* is an afhne parameter of the null 
geodesic vector defined by 



e = du. 

We complete the coordinate system by 



(B2) 

(B3) 

(B4) 

(B5) 
(B6) 



The * null tetrad related to this coordinate system can be 
given in terms of the original tetrad by a transformation 
of the form 



* 

u 


= u, 




r 


r 


= ^+ro 


c 


= c, 


c 


= c 



t = £, 

m* — 7 (m + Tl) , 
n* = n + Tfh + Tm + TTi. 



One can show that in our case 7=1 and 

r = -Igv-A. 



(B7) 
(B8) 
(B9) 



(BIO) 



The null tetrad {(* , n* ,m* ,in*) have the characteristic 
that is associated to a system {u* ,r* ,C,*), were r* an 
affine radial coordinate. Then we can use the known 
behavior of the spin coefficients in this tetrad-coordinate 
system [^^l in order to relate the leading order behavior 
of this quantities with those obtained in the more general 
tetrad (Z,n,m,m). 

In order to compute ctq, we must know the leading 
order behavior of t,(3, e', p' and n' . 

The spin coefficients transform according to 



which implies 



a = CF, 

T* ^T + Tp + fa; 

T = T* - Vp* ~fa*. 



(Bll) 
(B12) 
(B13) 

(B14) 



From this one can deduce that the leading order behavior 
of the original r has a 1/r term, namely 



1 



Similarly one has 

/3* = /? + f cr + re; 

with 

but our choice of tetrad is such that 

e = 0. 

Then in leading order the betas coincide, that is 



But one has 
and using that 
one obtains 



V* dP^ 1 

'71 ac^ ~ 2V" 



V 



00 = 



/3ov , 1 



with 



A 



A 2A^ 
V dPo 1 



V2 dC 2V 



Then one has 



Po = Pov + ^^vA. 



(B16) 
(B17) 
(B18) 

(B19) 

(B20) 

(B21) 
(B22) 
(B23) 
(B24) 



T-l 



dvA. 



(B15) 



Let us also note that 
- e'* = -e' - r/3' + f (r + /?) + T^a + Pf p; (B25) 

and 

- p'* = -f3' + Tp. (B26) 

Using that in the first two orders in the expansion of e' 
one has 

one can deduce that 

nine) + 2nCe = fln^iC); (B28) 

which was used above. 

Finally one can see that the other spin coefficients 
needed in order to compute a have the following behavior 

p'* =p + 2rf e - Tt' + f ^ct - f (/?' - (3) 
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and which imphes equation ((68)) . 

k' = nm{U); (B30) 

then one can see that they do not contribute to ctq. 

Collecting all these quantities in the expression for a, 
one has to leading order 

(B31) 



ao _ {wBf (l-^Y - 
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